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Abstract 

The multiple hadron production in the events induced by the heavy 
primary quarks in e~^e~ annihilation is reconsidered with account of 
corrected experimental data. New value for the multiplicity in bb 
events is presented on the basis of pQCD estimates. 



1 Introduction 

The so-called "naive model" [H [2] was the first attempt to give a framework 
for calculating the multiplicity of hadrons produced in addition to decay 
products of the heavy quark-antiquark pair in e~^e~ annihilation. Later on, 
it was argued [3] that the difference between multiplicities in heavy and light 
quark events (/ = u, d, s), 

5qi = Nqq{W)-Nu{W), (1) 

tends to a constant value at high collision energy: 

Sqi ^ = 2nQ - Nuimle). (2) 
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Here and in what follows, Nqq and A*";;- are mean multiplicities of charged 
hadrons in heavy and light quark events, respectively0 

The comparison with the data has shown that the "naive model" describes 
the data on 6bi up to = 58 GeV [H IH 13 [6] , but underestimates the LEP and 
SLAG data [II El IS]. As for the so-called MLLA formula ([2]), it significantly 
overestimates both low and high-energy data on 6bi- 

The detailed QCD calculations of the difference between associated mul- 
tiplicities of charged hadron in e~^e~ annihilation were made in [10]. The 
QCD expressions for 6qi from Ref. |T0] appeared to be in a good agreement 
with experimental measurements of associated hadron multiplicities in e+e~ 
annihilation (see, for instance, [II1[I2|)- Note that up to now, our formula 
provided the best description of all the available data on 6bi, see Fig. [TJ 

Moreover, we made a prediction for Sd [101 • It is also in a very good 
agreement with all the data on Sd [H El E] . Let us stress that the very value 
of 6d was derived in [10] before the precise measurements of 5d were made [8], 
that allows to test QCD calculations^] 

As we will see below, it is the hadron multiplicity in the light quark 
events that enables one to calculate the multiplicity differences 5qi. The 
mean charged multiplicities in // events at different energies corrected for 
detector effects as well as for initial state radiation were recently cited in [13] . 
The corrected multiplicity differences averaged over all presently published 
results were also presented 



(5^7^ = 3.12 ±0.14, (3) 
C'' = 1.0 ±0.4. (4) 

The first goal of this paper is to re-estimate our QCD predictions for the 
quantity 5^, taking into account the corrected experimental data on NijiW) 
from . The second goal is to argue that the MLLA formula ([2]) is nothing 
but some part of our QCD expression (see Section [2]), and, as the comparison 
to the data shows, it should be regarded as a rather rough approximation of 
the QCD result I 



^Everywhere below, it is assumed that we deal with mean multiplicities of charged 
hadrons. 

^In the low energy measurements [HIS], the total error of 5ci was about ±1.5. 
■^The shortcomings of the MLLA formula were already briefly discussed in Ref. [10] . 
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Figure 1: QCD prediction [10] and MLLA result [3] vs. experimental data 
on the multiplicity difference 5m- The data are not corrected as in Ref. [13]. 
The prediction of the "naive model" is also shown. 

2 QCD formula for multiplicity difference 

The hadron multiplicity in a qq event, Nqq{W), looks like [TO] 

Y 

Ngg{Y) = 2nq + Jdr]ng{Y -r])Eq{r]) , (5) 







where variables 



= In — (6) 
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and 

are introduced. In what follows, the notation q = Q will mean charm or 
beauty quarks, while the notation q = I will correspond to a massless case 
(when a pair of u, d or s-quarks is produced, whose masses are assumed to 
be zero). 

The first term in the r.h.s. of Eq. ([5]), 2nq, is the multiplicity of primary 
hadron decay products. It is extracted from the data (2nc = 5.2, 2ni, = 
11.0 [3], and 2ni = 2.4 pi]). 

The term Eglk"^ /W"^) in ([5]) is the inclusive spectrum of a gluon jet with 
a virtuality up to k"^ emitted by primary quarks0 It is defined by the dis- 
continuity of of the two-gluon irreducible 'j*g* {Z*g*) amplitude normalized 
to the total e"'"e~ rate. The quantity fig^k"^) is related to Ugi^k"^), the mean 
multiphcity of hadrons inside this jet: 

ng{k^) = ^^^n,ie). (8) 

Here as{k^) is a strong coupling constant, and Cp = {N^ — l)/2Nc, with Nc 
being the number of colors. 

The physical meaning of the function 

Y 

N,iY) = JdvhgiY-v)E,iv) (9) 



in Eq. ([5]) is the following. It describes the average number of hadrons pro- 
duced in virtual gluon jets emitted by primary quark (antiquark) of the type 
q. In other words, it is the multiplicity in qq event except for the multiplicity 
of the decay products of these quarks at the final stage of hadronization (the 
first term in ([5])). 

For the massless case, the function E = Ei was calculated in our pa- 
per [TU]. In terms of a dimensionless variable 

a = exp(— r/) (10) 



^It was explained in detail in Ref. [lOj that one should not consider this mechanism of 
hadron production via gluon jets as due to "a single cascading gluon" , as some authors 
believe [T^. That Eg is an inclusive spectrum of the gluon jets is seen, e.g., from the fact 



that J{dkyk^)Eq{kyW^) > 1. 
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it looks like 



E{7]{a)) = {l + 2a + 2a')ln--^-^{l-a)-a{l + a)(ln-] 

a 2 \ cr / 

+ 4a(l + (t)/((t) , (11) 

with 

1 

I{a)= [-^ln-^^-U,{l+a), (12) 
J I + X X 4 

(T 

where Li2(z) is the Euler dilogarithm. The function E{ri) is presented in 
Fig.d It has the asymptotics E{r])\^^oo = E^^'^^^r]) = r] - 1/2. 




Figure 2: The function E{y = t]). 

The derivative of E{ri) is positive, as one can see in the next Fig. [31 with 
dE{vi)/drj = at 7] = 0, and dE{rj)/drj = 1 at = oo. As a result, associative 
multiplicity Nq{W) ([9]) is a monotone increasing function of the energy W 
for any positive function ng(A;^)lfl 

Now let us consider the difference between multiplicities in heavy and light 
quark events, 5qi, which is defined by Eq. ([1]). The following representation 
was found in Ref. |lUj : 

^gf"^ = 2{nQ - m) - AiVQ(F„) . (13) 

^It results from the relation dNq{Y)ldY = J^drjUgirj) dE{Y - r])/dY. 
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Figure 3: The function dyE{y = rf). 



Here new notation, 



ANQ{Ym)=Ng-NQ= j dyn,{Y„,-y)AEQ{y) , 



as well as variables 



and 



F = In — 



(14) 

(15) 
(16) 



are introduced. 

The non-trivial result which was obtained in Ref. \XQ\ is that the function 
AEq = El — Eq = E — Eq (17) 
depends only on a variable 

p = exp{-y) , (18) 

but not on energy W. The explicit form of AEq is known to be 

7 1 20 

AEgiyip)) = (l-3p+-/)ln- + p(7p-20)J(p) + — -[1- J(p)] 

z p P ^ 



. 9 
+ 7p + - 



(19) 



where 



p-4 



J{P) 



-r— arctan , 



1=^ 



p > 4 
p = 4 
p < 4 . 



(20) 



Since AEq^i/) has the asymptotics 



- ^ exp(-|?/|) 



(21) 



the integral in Eq. (I14p converges rapidly at y — oo. The function AEg^y) 
is shown in Fig. H We find that AEq(?/)|j^^oo = AEl^"^"'\y) = y- 3/2. 

AEQ(y) 




Figure 4: The function AEg^y). 

One should mention the following important relationl^ 
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AEgiy - 1) - Eiy) 
In other words. 



^/e\n2 exp(-y/2) 



AEQ(y)^E(y + i; 



^Here (and below) e means the base of the natural logarithm. 



(22) 
(23) 
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If one puts AEgdi) = E{y + 1), then (neglecting also the contribution 
from the region y < —1): 



ANQ = Na{mle)-2ni. 



(24) 



Correspondingly, the approximate expression for 6qi is then of the form: 



Aappr) 



2nQ - Na{mQe) 



jrMLLA 



(25) 



where Sqi^^"^ is the MLLA prediction for the multiplicity difference from 
Ref [3]. Remember that the function Nii{W) describes the hadron multiplic- 
ity in light quark event at colliding energy W. 

However, Eq. (1231) is very far from being satisfied at relevant y < 
as it is clearly seen in Fig. [51 As a result, there could be a large difference 
between 6^J^"^ ([25]) and QCD expression ([13]). 



AEQ(y-l)-E(y) 
1 




8 10 



Figure 5: The difference AEqly — 1) — E{y) as a function of the variable y. 
To demonstrate this, it is convenient to represent expression for ANq ([T^ 



^For tlie beauty case, one has Ym < 3.2. 
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in the form: 



Ym + l 



/^NQ{Yra) = j dyhg{Ym + l-y)E{y) 



-1 

+ j dyUgiYm- y) AEgiy) 



-oo 



J dyngiYrr, + l-y)[AEQ{y-l)~Eiy)] 



^ [%(m|e) -2ni] + 5iV« {YJ + 6N^^^ {YJ , (26) 
that results in the formula (see Eq. (fT3|) ) 

5gf^ = 2nQ-Nuimle)-6N^^\Ym)-6Ni^\Y^) 

= 6^^r^-6N^^\Yj-6N^^\Yj. (27) 
Here we have introduced the notations 



-1 



N^^HyJ = J dyhg{Y^ - y) AEgiy) , (28) 



and 

Ym + l 



Nl^\ym)= J dyhg{Y^ + l-y)[AEQiy-l)-Eiy)]. (29) 



Note, both Nq (Ym) and Nq (Ym) are positive functions, since AEg^y) > 
at all y and AEqi^y - 1) - E{y) > at ?/ ^ (see Fig. 3 and Fig. 4). 
In order to exploit the corrected data on NijiW) at = 8 GeV, 

^"^^-(8.0 GeV) = 6.70 ± 0.34 , (30) 

we have chosen the mass of b-quark to be rrij, = 4.85 GeV, which corresponds 
to rrihy/e = 8 GeV. 

The estimates show that the dominant correction to Sqf^ is SNq \ not 
SNq^ . To calculate a lower bound of SN^"^^ , let us use the following inequality: 

AEgiy) = E{y + Ayg) . (31) 
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Note that Ayg is a monotone non-increasing function of ?/ ^ and it tends 
to 1 at large y. It solves the equation: 

^Eq{Y^) ^ E{Y„, + Ayg) , (32) 

where Ym is defined above f|T6|) . Then we get from Eqs. fl29l) and fl3Tl) : 

^iVg^ ^%(i;„ + AyQ)-%(r^ + l)- y rfyr2,(r„ + AyQ-y)i?(y). (33) 



For our further estimates, we need to know the hadron multiplicity in 
light quark events in the energy interval 2.5 GeV ^ ^ 28 GeV. By fitting 
the data on hadron multiplicity in the light quark events at low IV, we get 
the expression: 

NiiiW) = 2.07 + 1.11 \nW + 0.54 In^ W . (34) 

Putting Qq = 1 GeV, we find Ayb = 1.61. Taking into account that the 
last term in Eq. (!33l) is negligible jf| we get from (133|) . (!34l) : 

5A/'f ^ ^ 1.07 . (35) 

Correspondingly, our prediction accounting the revision of the data on the 
multiplicity in the light quark events, 

^ 2nb - NaiXm + Ayb) = 3.33 ± 0.38 , (36) 

appears to be lower than our previous result 6u = 3.68 [lOj. We used the 
value 

2^6 = 11.10 ±0.18 . (37) 

The error of A*";;- was taken to be ± 0.34. Let us stress that our upper bound 
is very close to the present experimental value of fl51) . 
Now let us derive a lower bound on 5^'"^. To do this, let us start from 
Eq. (El). It is convenient to represent the integral in (fT^ as a sum of two 
termso 

-1 n 
ANb = j dyhg{Yb-y)AEb{y) + j dy hg{Y„, - y) AE^iy) 

-4 -1 

= AN'^^ + AiVf ^ , (38) 



^Since E{y) < 0.02 in the region y ^ Ay,, - 1 0.61. 

^We took into account that the region — oo < y < —4 gives a neghgible contribution to 
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with Yfe = \n{ml/Ql) — 3-16. Consider the first term in (1551) . One can check 
that 

AE{y) < 0.18 E{y + 5.8) (39) 
in the region —4 < y < —1, that leads to the inequality 

4.8 

AiV(^) < 0.18^ dyng{Y, + 5.8 - y) AE,{y) . (40) 

1.8 

The estimations show that UgiYb + 5.8 — y) < 2 ^£,(4.8 — y) when y varies 
from 1.8 to 4.8. Thus, we get: 

AiV^^^^ < 0.36 [NiiiW = 11 GeV) - NaiW = 2.5 GeV)] = 1.54±0.17 . (41) 

The second term in (138!) can be estimated by using the inequality 

AE(y) < 0.62 E{y + 3.5) (42) 

which is valid in the region —l<y<Yb. Then 

AA/'f ^ < 0.62 [Na{W = 28 GeV) - NuiW = 3.5 GeV)] = 4.61 ± 0.30 . (43) 

As a result, we obtain from Eqs. (IT^ . (HM and (HTj) . the lower bound 

xQCD 
on (5^ : 

S^i^'^ > 2.55 ± 0.39 . (44) 

Fig. [6] demonstrates that our QCD predictions for 6bi are very close to the 
corrected experimental data. 

Our results can be compared with the MLLA expectation reported re- 
cently in Ref. [13]: 

^MLLA ^ 4 4 ^ _ (^45^) 

Note that the scheme of Ref. |13] is not stable against next-to-MLLA cor- 
rections. According to Eq. A(30) from [13], the MLLA prediction (123]) is 
modified as follows: 



SasiniQ) 



/TT^ 5 



5--- = 2n^ - N,imt,e){l + ^^^^ + ly " ijj | ' (^6) 
The next-to-MLLA corrections in (H6|) change the result (H5|) to 

^A^MLLA ^ 2.6 ± 0.4 . (47) 
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Figure 6: Our new QCD result for 6m (a corridor between the two solid lines) 
as well as MLLA prediction from Ref. [13] (a corridor between the two dashed 
lines) vs. corrected experimental data. 

The situation is worse in the case of c-quark. The formula fH6l) results in a 
unsatisfactory low value 

5™^^ = -0.1 ±0.4. (48) 

This demonstrates us once more that the lowest-order MLLA expression ([2]) 
is not correct. 

Moreover, as we have shown above, the deviation of the function AEqIu) 
from AEQ^^^\y) = y—3/2, as well as the deviation of E{y) from £'('^*f'")(?/) = 
y — 1/2, cannot be neglected. In other words, the MLLA formula (p5l) is, in 
fact, not a full QCD results It is nothing but a part of the correct QCD 
formulae ([T]), (HM in a very rough approximation E{y) = AEq^y — 1). So it 
is senseless to try to "improve" it with next-to-MLLA calculations. 

^°It explains why formula (|45p overestimates the data by more than one unit. 
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3 Conclusions 



We have derived the QCD formula for the difference between hadron multi- 
phcities in heavy and hght quark events in e^e~ annihilation (with Q being 
a type of a heavy quark): 



2?^Q - Niiimle) 



(49) 



<'0 
oo 



m: 



Q 



— n,(A:^)AEQU^ 



E 



fc2 



Here ng{k'^) describes the mean number of charged hadrons in the gluon jet 
with the virtuality up to fc^, and E, AEq are known functions. 

By using the data on the hadron multiplicity in light quark events Nii, 
corrected for the detector effects and initial state radiation effects ^3 
have obtained from (149!) the bounds: 



we 



2.2 < < 3.7 . 



(50) 



Let us note that this estimate does not depend on a specific choice of the 
function Ugik"^), and it is in a good agreement with the average experimental 



value = 3.12 ±0.14 



Two last terms in fH9|) are positive and numerically large0 As a result, 
a deviation of the MLLA prediction, 



^MLLA 

hi 



2nb - Nii{mle) 



(51) 



from the QCD expression)^ 



cQCD 



2K -^i)-J^ ng{e) ^Eq , (52) 

Ql 



^^In particular, the second term in (|49)) (dominating the third one) is equal to 1.1 for 
the case of the beauty pair production (toq = m;,, uq — ni,). 

^^This formula is an equivalent compact form of Eq. (f49|) for Q = b. 
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appears to be significant. 

As one can see, the MLLA formula is a too rough approximation of the 
QCD formula. The former results from the latter on the assumption that the 
quantities E{y) and AEq^y) can be replaced by their asymptotics ai y —>■ oo. 
Since the relevant values of y are far from being very large, this assumption 
is not correct, and it leads to a significant overestimation of Su- Thus, any 
attempt to use the MLLA expression ([2]) as a first-order approximation for 
higher-order calculations (as it is done in |13j ) is poorly justified. 
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